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Abstract 

Multi-species reaction-diffusion systems, with more-than-two-site inter- 
action on a one-dimensional lattice are considered. Necessary and suf- 
ficient constraints on the interaction rates are obtained, that guarantee 
the closedness of the time evolution equation for E%(t)'B, the expectation 
value of the product of certain linear combination of the number operators 
on n consecutive sites at time t. 

PACS numbers: 05.40.-a, 02.50.Ga 
Keywords: reaction-diffusion, generalized empty-interval method, 
next-nearest-neighbor interaction, multi-species 



1 mohamadi@azzahra.ac.ir 
2 mamwad@mailaps.org 



1 Introduction 



The study of the reaction-diffusion systems, has been an attractive area. A 
reaction-diffusion system consists of a collection of particles (of one or several 
species) moving and interacting with each other with specific probabilities (or 
rates in the case of continuous time variable). In the so called exclusion pro- 
cesses, any site of the lattice the particles move on, is either vacant or occupied 
by one particle. The aim of studying such systems, is of course to calculate the 
time evolution of such systems. But to find the complete time evolution of a 
reaction-diffusion system, is generally a very difficult (if not impossible) task. 

Reaction-diffusion systems have been studied using various methods: ana- 
lytical techniques, approximation methods, and simulation. The success of the 
approximation methods, may be different in different dimensions, as for exam- 
ple the mean field techniques, working good for high dimensions, generally do 
not give correct results for low dimensional systems. A large fraction of analyt- 
ical studies, belong to low-dimensional (specially one-dimensional) systems, as 
solving low-dimensional systems should in principle be easier [1-14]. 

Various classes of reaction-diffusion systems are called exactly-solvable, in 
different senses. In [15-17], integrability means that the A^-particle conditional 
probabilities' S-matrix is factorized into a product of 2-particle S-matrices. This 
is related to the fact that for systems solvable in this sense, there are a large 
number of conserved quantities. In [18-27], solvability means closedness of the 
evolution equation of the empty intervals (or their generalization) . 

The empty interval method (EIM) has been used to analyze the one dimen- 
sional dynamics of diffusion-limited coalescence [18-21]. Using this method, the 
probability that n consecutive sites are empty has been calculated. For the cases 
of finite reaction-rates, some approximate solutions have been obtained. EIM 
has been also generalized to study the kinetics of the g-state one-dimensional 
Potts model in the zero-temperature limit [22]. 

In [23] , all the one dimensional reaction-diffusion models with nearest neigh- 
bor interactions which can be exactly solved by EIM have been studied. EIM 
has also been used to study a model with next nearest neighbor interaction [24] . 
In [25], exactly solvable models through the empty-interval method, for more- 
than-two-site interactions were studied. There, conditions were obtained which 
are sufficient for a one-species system to be solvable through the EIM. In [26], 
the conventional EIM has been extended to a more generalized form. Using 
this extended version, a model has been studied, which can not be solved by 
conventional EIM. 

In a recent article [27] , we considered nearest- neighbor multi-species models 
on a one-dimensional lattice, we obtained necessary and sufficient conditions 
on the reaction rates, so that the time evolution equation for E* n {t) is closed. 
Here E* n (t) is the expectation of the product of a specific linear combina- 
tion of the number operators (corresponding to different species) at n consec- 
utive sites beginning from the fc-th site. All single-species left-right symmet- 
ric reaction-diffusion systems solvable through the generalized empty-interval 
method (GEIM), were classified. In this article, the method introduced in [27] 
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is generalized to more than two-site interactions. 

In section 2, multi-species systems with three-site interactions are investi- 
gated, which are solvable through the GEIM. Conditions necessary and suffi- 
cient for closedness of the evolution equation of E* (t) are obtained, and the 
evolution equation is explicitly given. 

In section 3, the more general case of multi-species systems with fc-site inter- 
actions, solvable through the GEIM, is investigated. Again, conditions necessary 
and sufficient for closedness of the evolution equation of E% (t) as well as the 
evolution equation itself, are obtained. 

Finally, in section 4, as an example a specific model with a three-site inter- 
action is introduced and exactly solved. 

2 Three-site interactions 

Consider a periodic lattice with L+l sites. Each site is either empty, or occupied 
with a particle of one of p possible species. Denote by N", the number operator 
of the particles of type a at the site i. a = p + 1 is regarded as a vacancy. N" 
is equal to one, if the site i is occupied by a particle of type a. Otherwise, iV" 
is zero. We also have a constraint 

s a N? = 1, (1) 

where s is a covector the components of which (s Q 's) are all equal to one. The 
constraint Q), simply says that every site, either is occupied by a particle of 
one type, or is empty. A representation for these observables is 

ATf := 1® — ®1 ®AT® 1® - - -(gil , (2) 

i-1 L+l-i 

where N a is a diagonal (p + 1) x (p + 1) matrix the only nonzero element of 
which is the a'th diagonal element, and the operators 1 in the above expression 
are also (p + 1) x (p + 1) matrices. It is seen that the constraint Q can be 
written as 

b-N=1, (3) 

where N is a vector the components of which are N°"s. The state of the system 
is characterized by a vector 

p g y ® ■ • ■ ® v , (4) 

L + l 

where V is a (p + l)-dimensional vector space. All the elements of the vector P 
are nonnegative, and 

S-P = l. (5) 

Here S is the tensor-product of L + 1 covectors s. 
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As the number operators N" are zero or one (and hence idempotent), the 
most general observable of such a system is the product of some of these number 
operators, or a sum of such terms. 

The evolution of the state of the system is given by 

P = HP, (6) 

where the Hamiltonian TL is stochastic, by which it is meant that its non-diagonal 
elements are nonnegative and 

SH = 0. (7) 
The interaction is a next-nearest-neighbor interaction: 

L+l 

H = J]] -ffi,i+i,i+2, (8) 

i=l 

where 

#M+M+2 := 1 ® • • • ® 1 ®H <g> 1 <g> • • -®1 . (9) 

z-1 L-i-l 

(It has been assumed that the sites of the system are identical, that is, the system 
is translation-invariant. Otherwise H in the right-hand side of would depend 
on i.) The three-site Hamiltonian H is stochastic, that is, its non-diagonal 
elements are nonnegative, and the sum of the elements of each of its columns 
vanishes: 

(s®s®s)# = 0. (10) 
Now consider a certain class of observables, namely 

k+n-l 

£l n := [] ( a ' N *), (11) 

l=k 

where a is a specific (p+ l)-dimcnsional covector, and Nj is a vector the compo- 
nents of which are the operators N". We want to find criteria for H, so that the 
evolutions of the expectations of £* n 's are closed, that is, the time-derivative of 
their expectation is expressible in terms of the expectations of „'s themselves. 
Denoting the expectations of these observables by n , 

El n := S £l n P. (12) 
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For 1 < n < L, we have 

n 

i=-i 

+ S £fe,n#fc-l,fc,fc+l P 
ri-2 

+ X! ^ £k,n H k-l+l,k+l,k+l+l P 

i=i 

+ S £k. n Hk+n-2,k+n-l,k+n P 

+ S £^ n Hk+ n -l,k+n,k+n+l P- (13) 

In the right-hand side of the last equality, the first two and the last two terms 
are contributions of the interactions at the boundaries of the block, while the 
summation term is the contribution of the sites in the bulk. 
To proceed, one may use the following identity 

s(b-N)=b, (14) 

using which, one arrives at 

(s (g s <g> s) [(a • N) ® (a • N) <g (a • N)ii] = (a ® a ® a)H. (15) 

Demanding that the bulk terms in (|13fl be expressible in terms of the expec- 
tations of n 's themselves, one arrives at the following condition for the bulk 
terms. 

(a<ga<ga)il = A(a (g> a (g> a), (16) 

for some constant A. By similar arguments, the condition coming from the left 
boundary terms of (|13(l can be obtained: 

(s <g s ® a (g> a)[(l <g> H) + (H ® 1)] = 

/i L s <g s (g s (g s + f L s (gi s ® s (gi a + 7 L s (g s (g a <g> a 

+ (5 L s <g> a (g a (g a + </> L a (g a <g> a (g a, 

(17) 

and finally, the condition coming from the right boundary terms of <|13[) is 

(a (g a <g> s (g s)[(l <g iJ) + (i? (g 1)] = 

/iRS (g s (g s <g s + v^a (g s (g s (g s + 7R,a (g a <g s (g s 

+ S^.a (g a (g a (g s + <g a (g a (g a, 

(18) 
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where the multipliers in the right-hand sides of (|17fl and l|18fl are constants. 
Arranging all these together, the evolution equation for E* n , for 1 < n < L, 
becomes 



E k,n =^ E k+2,n-2 + MR^fc,n-2 + v L E k+l.n-l + ^R^fc,n-1 

+ [tl + (n - 2) A + Tr]^ 

(19) 



In general, even if the initial conditions are not translationally invariant, one 
can solve the above equation. However, assuming that the initial conditions are 
translationally invariant, simplifies the calculations. Assuming that the initial 
conditions are translationally invariant, as the dynamics is also translationally 
invariant, E% n is independent of fc, and so one arrives at 



K = (ml + Mr)££- 2 + K + vr)K-i + + (« - 2)A + 7R ]£; 



The cases of zero-, one-, L-, and (L + l)-point functions should be considered 
separately. The zero-point function is equal to one. In the case of the one-point 
function, one arrives at an additional condition 



(s ® s <g> a ® s (g> s)[(l (g> 1 <g> H) + (1 g) H <g> 1) + (H <g> 1 <g> 1)] = 



Mo (s g> s g) s (g s <g s) + /i4(s®s®a®s®s 
+Ms(s (8>a(g)S(g)S®s) + Mi6(a ®s<g>s<g>s(g)s^ 
+Mi 2 (s <8>a(8ia(8)S(8)s) + M24(a (g>a®s®s(g)s^ 
+/X2s(a (g) a g) a g) s (g) s) + M2(s <g s (g) s <g a <g 
+M6(s <gsg)agiag)s) + Mi4(s g>a<ga<ga<gs^ 
+/i 7 (s 8s0a®a0a) + M3O 3 (gsgisgiatga^ 



In the case n = L, one arrives at the following additional condition 

(a <g> a® s <g> a ® a)[(l <g> 1 <g> if) + (1 <g> H <g> 1) + (i? ® 1 <g> 1)] = 



i/ 3 i (a ® a ® a ® a ® a) + i/ 2 7 (a ® a <g> s ® a ® a' 
+j/ 2 3(a (gs<ga<gagia) + fis(s ®a®a(g>a(g)a^ 
+^19 (a (g>s(gsg)ag)a)-|-z/7(s(gsg)ag)a(gai 
+^3 (s (g s g) s (g a (g a) + 1/29 (agiagiagisiga 
+^25(a g)ag)sg>sgia) + ^i7(a gisigsgisiga^ 
+^24(a ®a®s®s<g)s) + ^ 2 s( a <ga<ga<gsgis^ 



+ («J L + £ R )K+i + (0l + <fri)££ 



(20) 




(22) 



(21) 



+j/3o(a (g a (g a (g a g) s). 



(23) 
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Then the evolution equation for is 

E l =(^3 + vn + v 2 a)EI_ 2 + (v 7 + u 19 + u 25 + is 2 $)El_ 1 

+ + viz + v 21 + v 2% + v m )El + t>3iEl +1 . (24) 

The evolution equation for the case n = L + 1 is simply 

El +1 =\(L + l)El +1 . (25) 

3 fc-site interactions 

Now, let's go further and consider the more general case of fc-site interactions. 
The problem of single-species fc-site interactions solvable in the framework of 
empty interval method was addressed in [25] . In that article necessary conditions 
for a system to be solvable through the EIM was obtained. Here we obtain 
the necessary and sufficient conditions for a system to be solvable through the 
GEIM. To fix notation, let's introduce 

A ijl : = s® l <g>a® J ' ®s®', 

B ijl : = a®' <S> s®- 7 <g> a®', (26) 

where 

b® 1 :=b<g>---(8>b. (27) 

V v ' 

i 

First consider the bulk terms in the time derivative of £*'s. Demanding that 
the bulk terms can be expressed in terms of the expectations of £ *'s, one arrives 
at the following condition for the bulk terms. 

A ak0 H = XA oka . (28) 

The condition coming from the left boundary of the block is 

^fe-i.fc-i,0[ 1 ®fc-2 q H + + .. . + H® l® fe " 2 ] = C LoA ml °, (29) 

m+i=2fc-2 

while the condition coming from the right boundary of the block 

A o,k-i,k-i ^®k-2 R + . . . + H 1 ® fc -2] = C§ m A olm . (30) 

771 ,i 

771 + l = 2k-2 

For the blocks of the length j with j < k — 1, one needs the condition 
^fc-i^ fc -ij 1 ^+fe-2 0jff + !>> + i?(g}1 ®j+fe-2j = ^ C mp lA mpl - ( 31 ) 

771 .p. i 
7T l + P + (=j+2fc-2 
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Finally for the blocks of the length (L + 1 — j) with j < k — 1, one needs the 
condition 

B*-W.fc-i[l®i+*-2g |fl - + ... + fl -g )1 (8i+*-2] = ^ D mpl B mpl . (32) 

m,p,I 

Now it is easy to write the evolution equation of E*. For fc — 1 < n< L + 2 — fc, 

2k-2 

Ef l = (n-k + l)\El + Y, Et- k+ i +n Ck-2-i,L0 

1=0 
2fc-2 



^ a -fc+i+riC(^i,2fc-2-(j & l<n<L + 2 — k. 

(33) 



i=0 



For < n < k — 1, 



n+2fc-2 n+2fc-2-p 

^» = H C OT) p,„ +2 fc-2-p-m^p, < n < fc - 1. (34) 

p— m— 

Finally, fori + 2 — fc<n<i + l, 

n+2fc-2 n+2fc-2-p 

£<n = ^ ^ f m,p,n+2S;-2-p-m-Ex+l-p) L + 2 — fc < n < L + 1. 

p=0 m— 

(35) 

One also has, 

K = 1, 

El +1 = {L + l)\El +1 . (36) 

Equations i|55 )l to P fr can be solved to obtain a certain set of correlation 
functions (-EyJ's). 



4 A three-site model with annihilation, as an 
example 

Here, a model with three-site interactions is considered. Each site may be 
occupied or be a vacancy. A particle at each site may be annihilated, and the 
annihilation rate depends on its neighboring sites. The interactions are 

with the rate Ai, 
with the rate A2, 
with the rate A 3 , 
with the rate A 4 . (37) 



AAA -> A9A 
AA% -> A%% 
%AA -> 00A 
0A0 -► 000 
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(A and denote an occupied site and an empty site, respectively.) Choosing 
a = (1 0), it is seen that this model fulfills the conditions needed for solvability. 
This means that the evolution equations for E n 's with E n — (Ni ■ ■ ■ N n ) are 
closed. Using the method introduced in previous sections, one can obtain E^s. 
The evolution equation for E n 's are 

E n =-[(n- 2)Ai + A 2 + A 3 )]E n + (A 2 + A 3 - 2A 1 )E n+1 , 2<n<L-l 
E x = -A 4 £i - (A a + A 3 - 2A A )E 2 + (A 2 + A 3 - Ai - A 4 )£ 3 , 
E L = -(A 2 + A 3 )E L + (A 2 + A 3 - 2A 1 )E L+1 

E L +x = -{L+^AxEi (38) 

To simplify the model, assume Ai =0. This means that a block of particles, 
can only be altered at the boundaries. Defining 

a : = A 2 + A 3 , 

: = A 4 , (39) 

H38|) is rewritten as 

E n = a(E n+1 - E n ), 2 < n < L 

E x = [3(E 2 - Ex) + {a- /3){E 3 - E 2 ), 

E L+1 = 0. (40) 

There are two cases. 

i) a = 0. In this case, all the E n 's for n ^ 1 are constant and 

Ex{t) = e-^[Ex(Q) - 2E 2 (0) + E 3 (0)] + 2E 2 (0) - E 3 (0). (41) 

ii) a ^ 0. In this case, -El+i is constant. Using this, one can obtain El, and 
then El-i, and E 2 - The result is 



E 



L-k 



'JL+1 



E 

3=0 



iL-k+j 



(aty 



(42) 



where 7fc's are arbitrary constants to be determined from the initial conditions. 
Ex is also easily determined from the second equation of (|40|) . 
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